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APPLICATION OP REMOTE SENSING 


FOR FISHERY RESOURCE 
ASSESSMENT AND MONITORING 


INTRODUCTION 

This is the twelfth of a monthly series of progress reports required by the 
statement of work for Skylab Project 240 entitled "Application of Remote 
Sensing for Oceanic Gamefish Assessment and Monitoring" under Contract 
No. T8217B. 

OVERALL STATUS 

Several remaining analyses are being completed and the final report is in 
preparation. 

RESULTS 

The Principal Investigator presented a paper entitled "Preliminary Results of 
Fisheries Investigation Associated with Skylab", at a remote sensing symposium 
at Willow Run, Michigan, on 15 April 1974. 

A plot was drawn of the dolphin catch and the water discontinuities observed 
in the aerial photography. This plot was similar in format to one made earlier 

'v 

of the white marlin catch relative to the water rips. .Neither plot substantiates 
(as far as white marlin and dolphin are concerned) an opinion held by fishermen 
that better fishing may be found in the vicinity of rips. 



Remotely inferred values for sea surface temperature, chlorophyll-*! and 
turbidity were substituted for sea truth measurements in prediction moders 
developed in previous analysis. Model performance, using the new values, was 
disappointing . 

EXPECTED ACCOMPLISHMENTS 

Correlations of dolphin with the environmental parameters will be explored in 
the time remaining on the project contract. First priority in analysis had 
been accorded to a billfish species. An analysis relative to white marlin is now 
complete. Dolphin were selected for follow-on analysis since more dolphin 
were caught than any other gamefish species during the tournament associated 
with the project mission phase. 

The draft final report will be. mailed during the coming month. May 1974, as 
required contractually. 
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ABSTRACT 


SCALE AND ORIENTATION CONTROL IN GEODETIC 
APPLICATIONS OF SATELLITE ALTIMETRY 


An analytical method for geodetic computation of the marine 
geoid (the geo id in the oceans) from satellite altimetry is developed and 
validated with data from Skylab mission SL-2. The criteria for achieving 
accurate scale and orientation of satellite altimetry geoid are shown to 
require marine geodetic control to offset systematic errors in the orbit 
(orientation is completely orbit dependent) and the altimeter data. 
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REPORTS AND DATA RECEIVED 
(Continued) 


(8) SKYLAB IV/EREP DATA BOOKS PROCESSED FROM 


DDC Accession No. 

34-15362 

34-15176 

34-05945 

34-15044 


Microfilm Roll No. 

34-07854 

34-17168 

34-17014 

34-17199 


DPAR No. 

S193B-69-6-54-72-2 
S193B-86 -2-85 -32-1 
S 190A- 106 -3-54-22-1 
S190A- 108- 1-97- 22-1 


(9) SKYLAB IV/EREP DATA BOOKS 
DDC Accession No. 


Microfilm Roll No. 


DPAR No. 


34-15408 

34- 

34-15407 

34- 

34-15412 

34- 

34-15411 

34- 


1 EARTH RESOURCES DATA 


17037 

S193B-69-7-85-73-B 

17055 

S193B-69-6-54-73-B 

17085 

S193B-71-1-97-72-A 

17172 

S193B-71-1-97-73-B 

FORMAL CONTROL 

BOOK 


No. of 
Copies 

4 


EREP Bulletin, Number 
From: FS5 6/Data 

EREP Bulletin, Number 
EREP Bulletin, Number 
EREP Bulletin, Number 


1 December 11, 1973 

Systems and Analysis Directorate, Data Manager 

2 December 11, 1973 (4 Enclosures) 

3 December 11, 1973 (with Cautionary Note) 

4 January 30, 1974 
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PROGRESS 


During this reporting period, our main effort involved 

(1) Completion of our investigation of scale and orientation 
control in geodetic applications of satellite altimetry. 
This is an outgrowth and a necessary integral part of 
our investigation of calibration and evaluation of 
Skylab altimetry for geodetic determination of the 
geoid, In particular, the accuracy problems of 

SKYBET necessitated this investigation. A formal 
write up of our results is given as Appendix B* 

(2) Continued work on the remaining SL-2 data 

(3) Continued investigation of effect of orbit errors 
in geoid determination from satellite altimetry 

(4) Quick look examination of SL-4 data tabulations 
received so far. 
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SCALE AND ORIENTATION CONTROL IN GEODETIC 
APPLICATIONS OF SATELLITE ALTIMETRY 


INTRODUCTION 


Determination of the geoid in the form of heights above a 
reference ellipsoid is the basic step in the geodetic applications of 
satellite altimetry. The geoid to be determined must be in absolute 
position or geocentric (i.e., centered at the earth's center of mass) and 
have correct scale, shape and orientation in order to meet the goals of 
geodesy and also make contributions to the solution of problems in earth 
gravity modeling, geophysics, oceanography, etc. Correctness of shape 
depends on the precision of the altimeter and, in theory, absolute centering 
and orientation are dependent on the satellite orbit ephemeris. The 
correctness of geoid scale requires that the orbit ephemeris and the 
altimeter either have no biases or systematic errors, or that such biases and 
systematic errors must be known to an accuracy better than the error tolerance 
of the geoid to be computed. Currently and for sometime to come, these two 
scalar conditions cannot be met because of unknown systematic errors or biases 
in tracking station geocentric coordinates, the earth's gravity model, the 
tracking systems and the altimeter itself. There is , therefore, a need for 
other sources of scale and orientation control. This paper discusses the 
need for and the use of terrestrial marine geodetic data to obtain scale and 
orientation control in the computation of the marine geoid (i.e., the geoid 
in the ocean areas) from satellite altimetry. 

Three types of terrestrial geodetic parameters are required for 
this scale and orientation control: (1) The best available estimates of the 

figure of the earth in terms of the size and shape of a reference ellipsoid, 
(2) geoid heights references to this ellipsoid, and (3) marine geodetic 
controls. The first two of these are required as a priori inputs to provide 
a coarse scale. The third serves as benchmarks establishing the fine scale 
and misclosure errors. This is akin to leveling practice on land. Satellite 
altimetry is simply geodetic leveling from space. 

Various estimates of the figure of the earth are in Mueller [1966] 
and Khan [1973]. The best space age estimates of the equatorial radius value 
range from 6,378,124 m. [Strange et al. 1971] to 6,378,169 m. [Veis, 1967], 
with most estimates in 6,378,140 + 5 m range and a flattening of l/298’255. 
Unfortunately, for the geoid, agreement and/or compatibility of various 
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authors geoids is considered very discouraging [Decker, 1973 and Fubara, et al* 
1972a], particularly in the ocean areas. Vincent and Marsh [1973] geoid based 
on equatorial radius of 6,378,142 m. and flattening of 1/298*255 was selected 
for this investigation. A comparison of marine geoids of Vincent and Marsh 
[1973], Vincent, Strange and Marsh [1972] (see Figure 3) and Talwani, et al. 
[1972] shows why the choice of any of them can provide only coarse scale. 

The fine scale requires the use of marine geodetic control established via 
the use of satellite geodesy and as trogravimetric techniques [Fubara and 
Mourad, 1972a, 1972b and Fubara, 1973a, Hour ad , et al. , 1972a and 1972b], 

Figure 1 shows schematic geocentric relations of the various surfaces 
associated with satellite altimetry. TM is the raw altimeter range which has 
to be corrected for laboratory instrumental calibration, electro-magnetic 
effects, sea state, and periodic sea surface influences to give TS. S 
represents the non-periodic n sea level". CT and CE, the geocentric radii of 
the altimeter and E, its subsatellite point on the reference ellipsoid, are 
computed from satellite tracking information. EG is the absolute geoidal 
undulation to be computed from this investigation, while SG is the quasi- 
stationary -departure of the mean instantaneous sea surface from the geoid - 
the "undisturbed 11 mean sea level. Details of the altimetry data from Skylab 
mission SL-2, EREP pass 9, used in this paper are in Fubara and Mourad [1974]. 


t Satellite altimeter 



FIGURE 1. GEOCENTRIC RELATIONS OF SURFACES 
INVOLVED IN SATELLITE ALTIMETRY 
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ANALYTICAL DATA HANDLING FORMULATIONS 


Each measured altimeter range R? with an associated measurement 

residual v. is intrinsically related to (1) X , V and Z (the satellite 
i s s s 

coordinates at the instant of measurement), (2) the absolute geoidal 

£ 

undulation N^ (of the subsatellite point) based on a reference ellipsoid of 
of parameters a (equatorial radius), and e (first eccentricity) and (3) the 
biases in all measurement systems involved. The condition equation for this 
intrinsic relationship can be stated as: 

v. + R°(l + Ac) - h. + N? + AN. = 0 (1) 

XI X 1 1 


where 


Ac = f. (systematic errors in X , Y , Z } the altimeter bias and 

X s* s s 

sea state correction bias) is the total geodetic calibration 
factor to be determined* This factor controls the scale or 
the radial relation to geocenter of the geoid deduced from 
satellite altimetry. 

a o o a 

N. = N. 4 AN. (N. is an approximate value for N, ) = f^(a, e): 
xx ii r * i 2 

and iu is the geodetic satellite height above the reference ellipsoid, or 

h i - f 2 <V V V 5 . 

where a and e are parameters of the reference ellipsoid for the geodetic 

datum of the tracking stations whose coordinates are used in computing the 

satellite coordinates X 3 T , and 2 * Equation (1) presumes that a - a 

s s s 

and e^e; and also th3t the two reference ellipsoids are concentric and 
geocentric . 

In current geodetic practice, because of multiplicity of geodetic 
datums and the non-existence of an universally accepted datum, the a =* a, etc. 
requirements are hardly ever met. A geodetic datum is uniquely determined by 
seven parameters. One such set of parameters is a, e, Ax, Ay, Ae , Ac, and ATI. 
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a and e define the size and shape of the reference ellipsoid; Ax, Ay and 

Az relate the center of the reference ellipsoid to the geocenter and are 

purely translatory; while A| and AT| are angular values to ensure parallelity 

between the minor and major axes of the reference ellipsoid and the mean 

rotational axis and mean terrestrial equator of the earth respectively. 

For each geodetic datum, every effort is made to ensure that Af = AT* = 0. 

However, as shown in Fubara and Mourad [1972a], this condition has never 

been exactly realized but its effect can be neglected. 

The change Ah in h due to the changes Aa and Af in the dimensions 
l i 

of the reference ellipsoid and Ax^, Ay^, and Az^ in its position relative 
to geocenter is given by Heiskanen and Moritz [1967] as 

2 

Ah =_ Cos’cCosX Ax - CoscSinXAy - SirgrAz - Aa + a Sin (pAf (2) 

i ' o o o 


where 


f = flattening of reference ellipsoid [f = 1 


(l-eV /2 } 


cp and X = geodetic latitude and longitude corresponding to X g , 

Y , and Z . 
s s 


To include the Ah^ correction parameter, Equation (1) should be rewritten as 


v + R (1 + Ac) - (h. + Ah.) + N° + AN. - 0 (3) 

i i i x x i 


to reflect changes in reference ellipsoidal parameters whenever necessary. 


Equation (1) or (3) can be rewritten in matrix form as 

F 1 X 2’ L l } = °* (4) 

subject to the normalized weighting functions and associated with 

and respectively. 

In this model, all parameters and measurements of the mathematical model are 
treated as "measurements 11 and weighted accordingly. This mathematical model 
for the generalized least squares processing of experimental data is based on 
works of Schmid and Schmid [1964]. 

XX* 
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The superscript n a u denotes the exact true values of the "'neasuremen ts ,f . 
Usually, these true values are not known. Instead, the corresponding 

measured or approximate values X° X°, and L with associated variance- 

_1 ^1 -1 i z i 

covariances P , F 9 , , are estimated or measured. Therefore, 

Equation (4) can be rewritten in the form 


F 2 [ (X 1 + A l ) ’ (X 2 + A 2 ) ’ (L 1 + V l>] = ° 

where in relation to Equations (1) and (3) 

x i = x i + A i =N ° + AN i 

X 2 = X 2 + A 2 = R i AC “ (h i + Ah i 


(5) 


The linearized form of Equation (8) is 

A A + B x A 2 + + F 2 (X°, X°, L°) = 0 (6) 

A , B , and C are the first partial derivatives in a Taylor series 
1 1 * o o o . 

expansion of Equation (5), associated with X^> X^, and L^, respectively, 

vhile A 2 > and V are the correction parameters to be determined. 

The least squares solution of Equation (6) to derive the corrections 

A , A , and V to "measured" X° X° and L° is as developed in Fubara [1973b]. 
12 1 1 z i 


DATA PROCESSING , RESULTS AND ANALYSIS 

The data being analyzed consist of two sets “ A and B. Set B 
altimeter ranges are the same as Set A plus some biases, and orbit B and 
orbit A differ by a set of systematic errors. 

The elementary deduction of geoid heights from satellite altimetry 
is to merely subtract each range measured by the altimeter from the corresponding 
geodetic height of the satellite computed from the orbit ephemeris. 

Compared to Equation (3), this simple form which is shown to have 
many deficiencies is 


R. - (h- + Ahi) + N i = 0 


(7) 
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Figure 2 gives the results of applying Equation (7) to the four combinations 
of data Sets A and B, and also a profile of the same segment from Vincent 
and Marsh [1973] geoid which is reliable to + 5 to + 15 meters in the oceans 
according to the authors* It is obvious, therefore, that as long as computed 
orbits and the altimeter data have systematic errors or biases, this simple 
subtraction technique cannot furnish the geoid with correct scale and 
orientation relative to geocenter. 

The results from processing the same data combinations, using the 
analytical formulations of Equations (3) to (6), and a priori geoid height 
input from Vincent and Marsh [1973] are shown in Figure 3. The match up 
between profiles AA and AB (same orbit but different altimeter data biases), 
and also between BB and BA (same orbit and different altimeter data biases) 
indicate (a) the analytical processing is efficient; (b) based on this type 
of processing, input of a priori geoid information does provide useful "coarse 
scale" to the satellite altimetry geoid irrespective of unknown altimeter 
system biases; and (c) marine geodetic controls (station whose geocentric 
three-dimensional coordinates are accurately known) are required to provide 
the necessary fine scale and absolute orientation unless the orbit is 
errorless and can be validated by other independent means. 

Orbit A and the Vincent and Marsh [1973] geoid were each computed 
from nearly identical sets of geopotential coefficient. Therefore, the 
excellent agreement between the satellite altimetry geoid segments AA and AB 
based on Orbit A and the Vincent and Marsh geoid is not a proof of the accuracy 
of either the orbit or that geoid. Rather it appears to indicate the precision 
of the Sky lab altimeter. 

The effects of errors in a priori geoid height inputs and the scale 
dependency of satellite altimetry geoid on geodetic control were investigated. 
The results are shown in Figure 4* Orbit A and Set A altimeter ranges were 
used and A-I, B-I and GG-73 were the a priori geoid height inputs. The 
corresponding "satellite altimetry geoid" profiles deduced were A-0, B-0 and 
AA. These results indicate that, provided the precision of the altimeter is 
reliable, errors in a priori geoid height input affect only the scale and not 
the shape of the resultant altimetry geoid. 
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Gcoid "Height " f meters 


Lat . 
Long. 
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36 . 90N 
286.04 


35.13N 

288.68 


35.73N 

291.20 


Latitude and Longitude in degrees 


■pTr*TTT5T7 O 

X i_ V U Iiu t 


CONVENTIONAL C-EOID AND SATELLITE -ALTIMETRY GEOIB 
SEGEMENTS (SKYLAB SL-2 EREP PASS 9 DATA) 





Geoid "Height", meters 


A-l * A Priori Geoid Input (Errors purposely added to GC-73) 
A-0 * Resultant Satellite Altimetry Geoid from A-I 
CG-73 = A Priori Geoid Ir.put (Ground Truth Data) 

AA = Resultant Satellite Altimetry Cecid from CC-73 
B-^I = A Priori Ceoid Input 

B-0 ~ Resultant Satellite Altimetry Geoid from B-I 


Latr 3G.90N* 
Long, 265. 04 e 


35. 1 3N 
. 268,6 Sr 

(Latitude (j Longitude in degrees) 


33.30N 
291, 20E 


FIGURE 4. EFFECT OF ERRORS IN A PRIORI GEOID HEIGHT INPUTS AND SCALE 
DEPENDENCY OF SATELLITE ALTIMETRY GEOID HEIGHT ON GEODETIC 
CONTROL (GROUND TRUTH) 
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CONCLUSIONS 

Based on these investigations, it can be said that: 

(1) Analytical data handling formulations of the type developed are 
effective and necessary for accommodating data biases to ensure 
a reliable scale in geodetic application of satellite altimetry* 

(2) A priori geoidal information is required to give a coarse scale* 

The use of geodetic controls or benchmarks, at which absolute 
geoid heights are known, is indispensable to ensure that the 
deduced geodetic heights are absolute, correctly scaled and 
oriented* The satellite orbit and the altimeter biases or 
systematic errors are known more accurately than the tolerable 
error of the computed geoid* The establishment of such controls 
by marine geodesy from a combination of as trogravimetry and 
satellite data is discussed in Mourad and Fubara [1972], and in 
Fubara and Mourad [1972a] and the practical implementation is 
partially demonstrated in Fubara and Mourad [1972b]* There is 

an implicit correlation between this conslusion and the conclusion 
based on a different type of investigation in Rapp [1971] that: 
n In carrying out simulation studies with non-global data it was 
concluded that altimetry data could not be used alone for potential 
coefficient determination.*** Consequently, the altimetry data 
was combined with geoid undulation information in non-ocean blocks 
and with existing terrestrial gravity data 11 . Therefore, marine 
geodesy appears indispensable for the full achievement of 
satellite altimetry objectives of GEOS-C, SEASAT and the NASA- 
proposed "Earth and Ocean Physics Applications Program" • 

(3) Accurate deduction of the geoid from satellite altimetry cannot 
be achieved by merely subtracting altimeter ranges from the 
corresponding geodetic heights of the satellite unless to present 
tolerance levels, (a) the satellite orbit is errorless, (b) the 
altimeter does not drift, and (c) the altimeter system biases 
are either non-existent or are absolutely known* Even then. 




the proof that preset tolerance levels are achieved will 
depend on marine geodetic control. Therefore, as of now, 
satellite altimetry ranges cannot be regarded as representing 
direct determination of absolute geoid heights* However, 
preliminary results from the Skylab altimeter show that 
satellite altimetry, backed by marine geodesy can accomplish 
EOPAP objectives. 
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